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ABSTRACT 


This report derives from elementary principles the general equations of 
motion for a missile utilizing a fixed-plane coordinate system, i.e., a co¬ 
ordinate system with one axis constrained to lie in a given plane. 

Included in the derivation are explicit expressions for introducing v»:nd 
and an alternate set of equations to be used when singularity conditions 
are approached. Means are provided for automatically converting to the 
alternate set of equations so that uninterrupted trajectory simulation can 
proceed under all conditions. A complete discussion of initial conditions 
is included. 

The general equations can be used for flat or spherical, rotating, or non¬ 
rotating earth cases. 


INTRODUCTION 


This report is part of a continuing program to give Picatinny Arsenal 
a complete capability in the flight simulation of all types of projectiles 
and missiles, whether ballistic or rocket-boosted, guided or unguided. 

Contained herein is a rederivation, from elementary principles, and an 
elaboration of the fixed plane coordinate system described in Reference 1, 
"Trajectory Equations for a Six-Degree-of-Freedom Missile." In Refer¬ 
ence 1, a list of only the basic equations unique to the fixed plane co- 
ordina e system is presented. This report completes this list and includes, 
for the first time, a derivation for initial conditions, as well as an alter¬ 
nate set of equations to be used when singularity conditions are approached. 
Means are provided for converting to this set during flight simulation. 
Finally, explicit equations for introducing wind into the equations of 
motion are derived. 

This alternate set of fixed-plane trajectory equations was required to 
supplement (and, in many cases, replace) the existing equations utilizing 
the missile-fixed coordinate system as derived in Reference 1. To be 
explicit, the previous equations produced satisfactory trajectories for low- 
spin projectiles, but were inadequate for the simulation of spin-stabilized 
shells. The most obvious differences were in the large deflections accom¬ 
panying most trajectories for high-spin projectiles, these being three or 
four times as large as had been anticipated. Discussion of this matter 
with personnel from the Naval Weapons Laboratory, Dahlgren, Virginia, 
indicated the existence of a narrow band of permissible integration step 
sizes (incremental time steps) whereby both truncation and round-off 
errors are of acceptable magnitudes. It is possible that when forces and 
moments are referred to a missile-fixed coordinate system (as was the case 
in Reference 1), and high-spin rates arc to be accounted for, this band of 
acceptable time increments becomes even more narrow or perhaps non¬ 
existent. Although this is not known with certainty, it provides the im¬ 
petus for studying the fixed-plane coordinate system described in this 
report. In particular, this coordinate system has one axis constrained to 
lie in a given plane and, consequently, does not rotate with the missile. 

To provide a working simulation and to determine whether this new co¬ 
ordinate system alleviates the conditions producing unsatisfactory trajec¬ 
tories, the equations derived in this report were programmed by the 
Digital Applications Unit for the IBM 709 computer, Reference 5 contains 
a description of the corresponding computer program. 
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The equations again consider both flat and spherical, rotating and non- i 

rotating earth cases, and again make use of the Euler transformations * ' 

(as opposed to direction cosines) to express vector quantities in various 
coordinate systems. 

As before, guidance factors, motion of the earth along its orbit, launcher 
effects, and asymmetric missiles are not included in the derivations. 

These represent potential areas of extension of the present equations of 
motion. 


NOTATION 

Unit vectors in each of the three orthogonal directions are represented 
by ij , j jt k; where subscript i denote^ the coordinate system under con¬ 
sideration. Components of vectors in each of these directions will have 
two subscripts. The first subscript (X, Y, or Z) denotes a component 
along the i, j, k axis, respectively; the second subscript (I, E, etc.) 
denotes the referencing coordinate system. Thus V XE is the component 

of the vector V along the i axis of the E-coordinate system. 

Arrows over vector quantities denote vectors of arbitrary magnitude, 
bars over vector quantities denote vectors of unit magnitude. 

Subscripted vectors other than the subscripts mentioned above are 
enclosed in parentheses, i.e. (<3^). 


Finally _ 


denotes time differentiation of vectors relative to the 


. th 
i 


dt 

coordinate system, while a dot over a given variable denotes scalar dif¬ 
ferentiation with respect to time. 
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SYMBOLS 


Vector Quantities 


‘i* >i* k i 
‘E’ >E* k E 

iw> k W 

: tf ! H’ k H 

*v> i V k v 
*M’ 'm 1 k M 

u 

» 

a 


oj 

Q 

(w T ) 

-e 

CJ* 

-* 

R 

e 

V 

(V w ) 

If 

li 

•e 

J 

9 


(V) 


coordinate systems 


Angular velocity of (H) relative to (I) 

Angular velocity of (E) relative to (I), i.e., angular velocity of 
the earth about its axis 

Angular velocity of (H) relative to (E) 

Angular velocity of (M) relative to (H) 

Angular velocity of (M) relative to (I) 

Angular velocity of (V) relative to (E) 

Vector from center of earth to the current missile CG position 

-a 

Time rate of change of R. i.e., missile velocity 
Vector describing wind 
Summation of forces acting on missile 
Summation of moments acting on missile 
Angular momentum of missile 

Vector giving direction of gravitational force exerted on missile 
Velocity of missile relative to the air 
Thrust mil alignment vector 


Thrust vector 
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SYMBOLS (Cent) 


Scalar Quantities 


o, o, e, v 

<jxj)' 

p 

d 

k r>A’ k N’ k F 
k H' k <^p’ k d> 

r c 

a n 

8 

R e»rth 

*o 

h 

*r *a 
k 

a, /I 

*v w l 

d> 

A, B, G, H 


mass of rocket 

longitudinal moment of inertia 
transverse moment of inertia 
Eulerian angles 

Angles relating (H) to (V) coordinates 
Density of air 
Diameter of missile 
Aerodynamic coefficients 

Aerodynamic coefficients 

Distance from nose of missile to CG 

Distance from nose of missile to normal center of pressure 
Distance from nose of missile to magnus center of pressure 

Gravitational acceleration 
Radius of earth 

Gravitational acceleration at sea level 

Altitude of rocket 

Thrust malalignment angles 

Portion of thrust devoted to produce jet tortjue 

Anglos defining missile position 

Magnitude of wind velocity 

Angle defining direction of wind 

Angles defining initial conditions 



SYMBOLS (Cont) 


. Denotes time differentiation 

" Denotes quantities referred to (V) coordinates 

>; Denotes the vector cross product 


PROCEDURE 


Several coordinate systems are utilized in deriving the equations of 
motion. They are tabulated here for later convenience. 

i|, jj, kj Inertial coordinate system; origin (O) at the center of the 
earth, k ( axis coincident with positive spin axis of the 

earth, ij axis coincident with 0° longitude, and jj axis so 
directed as to form a right-handed coordinate system. 

i, j e , k E Earth-Hxed coordinate sys em; identical position as 

ij, jj, kj at time equal zero, but coordinates are to be fixed 
to the earth. 


*W’ iw> Azimuthal coordinate system; this coordinate system is used 
to introduce wind data into the equations of motion. The 
origin of these coordinates is again at O, the i v axis is 

directed towards the missile CG, the k\y axis is directed 
towards the positive spin axis of the earth, and j w is so 
directed as to produce a right-handed coordinate system. 

*H’ Ih* Fixed horizontal plane coordinate system; i H axis along 

missile longitudinal axis is directed from the centroid 
towards the nose of the missile, is constrained to lie in 

a horizontal plane, parallel to the 7 F , plane, and k~ H is 

directed as to produce a right-handed coordinate system. 
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Fixed vertical plane coordinate system; the iy axis is coin- \| 

cident with the i^ axis, the j y axis is constrained to lie per- v 

pendicular to the i £ axis (i.e., remain in a specified vertical * 

plane), and the k v axis is so directed as to determine a right- 
handed coordinate system. 

Missile-fixed coordinate system; the iy axis is coincident 
with missile longitudinal axis, the and k M axes are rigidly 
attached to the missile to form a right-handed coordinate 
system. 

The procedure will be divided into several sections, as outlined in the 
Table of Contents, each describing a particular aspect of the equations of 
motion, with the last section combining all that precedes it and serving 
as a summary of the equations of motion. The first section will derive the 
general equations of motion, and will include the necessary transforma¬ 
tions between coordinate systems. 

The General Equations of Motion 

As usual, Equations 1 and 2, the basis for Newtonian mechanics, pro¬ 
vide the foundation for the equations of motion. 1 


*V > i V» k v 


‘M. Im. k M 


-♦ 

If 



dt J 


(i) 



dt 


( 2 ) 


Strictly speaking Kquation 1 should read: 


a 

V F 


A. I mV) 
dt 


d.V 

m _!_ * urn 


dt 


d’.R - . 

m _!_ i um 

dt*' 


<!') 


The latter term, un, involves the amount of mass being expelled from the rocket system, 
and the relative velocities of the rocket and the exist gases, (u>. This term is known as 
the jet reaction and is considered aspart of the thrust, whose total makeup also includes 

considerations of nozzle deaign, operating temperatures, and pressures. Kquation l‘ uses 
u in place of the seemingly indicated V becajse the rocket system itself hasnot really 
been defined in this presentation. The interested reader is referred to References 2 and 
} for a more complete discussion. 
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Here 2F and 2L are summations of forces and moments acting on the mis¬ 
sile, R is a vector from the center of the earth to the current CG (center of 

gravity) position of the missile, and J is the total angular momentum of the 
rocket. 

Since Equations 1 and 2 are vector equations in three-dimensional 
space, three component or scalar equations are implicit to each. To deter¬ 
mine explicit directions relative to the missile along which forces and 
moments can be conveniently summed, consider a coordinate system whose 
origin is at the missile CG and whose axes are directed as given previously 

by i, j k _ Let d be the angular velocity of this coordinate system 
H H n 

relative to inertial coordinates, where 


^ ' “xh‘h + “YH’H + U ZH k H ^ 

in terms of the (H) coordinates. 

Utilizing the derivation presented in Appendix A, Equation 1 may be 
rewritten as follows: 



where all quantities in the right hand side of the equation are understood 
to be expressed in (H) coordinates. One may correctly surmise that the 
presence of the « v term relates the motion of the coordinate system (H), 
to which forces and moments are referred, to an inertial coordinate system 
(I) as required by Newton’s Laws of Motion. 

Rather than differentiate Equation 4 directly, it is convenient to let 


* V XH i H < V YH j H t V 7H k H 
dt 

so that Equation 4 becomes 


It 



dt 


m 


(V 

\ dt 



(5) 


( 6 ) 
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Since V is already expressed in (H) coordinates, (see Equation 5), one 
can write 

d H V . - . - . — 

—— “ V XH*H + V YH>H f V ZH k H W 

dt 

Combining Equations 3, 5, 6, and 7, and performing the indicated opera¬ 
tions, one then obtains the three component force equations, namely: 

= m[V XH + (“ YH V ZH -«ZH V YH )1 in the r H direction 
= m[V YH +(W ZH V XH " W XH V ZH )] in th<? difeCtion (8) 

= + ^ W XH^YH ~ “ YH^XH^ in the k H direction 

One can treat the moment equation (Equation 2) in a similar manner, 
obtaining first 


IF 


XH 


IF 


YH 


IF 


ZH 


XL = ^ = iil r cS X f (9) 

dt dt 

Denoting (w x ) as the total angular velocity of the missile relative to 

inertial coordinates, we can write the angular momentum in the following 
form: 


J ^X (<J T } XH'h + I Y (w T ) YH i H f I Y (to T ) ZH k H 


( 10 ) 


where I x , I y are the moments of inertia of the missile about the longi¬ 
tudinal and transverse axes, respectively. Note that the assumption of 
rotational symmetry (I z = I y ) is implicitly made in Equation 10. 

Since J is already expressed in the (H) coordinate system, the time rate 
of J can be expressed as follows:' 


^X^ T^XH* H ' l y (w x )YniH 

dt 


I y (o. T )zH ,t H 


(ID 


Again, strictly speaking, terms of the folloatng form should be included in the differ* 
enti.ition during the burning stages: 

'^^Vxir *y (w t ) yh’ 'y^'t^zh 

Howevei, for a fairly stable rocket (no tumbling), l y ( T^YIf 10 ^ ^Y*°' , 1^ ZH c,n 
neglected. The expression Iy(o) X ) X j| can be written as follows: 

i x (w-^XH ink^or-^xH * m(0> T>XH <k,) 2 

where k is the radius of gyration. The latter term can be neglected because usually the 
burning fuel has Tittle effect on k\ In addition, for low-spinning rockets, the former term 
also can be neglected. For high-spin rockets, however, it k <"> x > x „ «* combined with 

other terms which comprise the jet torque which acts to produce high-spin rates. The 
reader is referred to Reference } for further details on this matter. 
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Appropriately combining Equations 3, 9, 10, and 11 now results in the 
component equations for the moments: 

^ L XH = ‘^T^XH + (“yhM^T^ZH ~ "ZhM^T^Yh) 

2L ™ = Iy(<y t) yh + ("zh^t) xh - "xhM^tW (12) 

SL ZH = M^ZH + <"xhM w ^YH " "YH^tW 

Equations 8 and 12 comprise the basic six-degree equations of motion. 

These equations are of limited use, however, until one knows explicitly 
how the (H) coordinate system is moving with respect to inertial or earth- 
fixed coordinates. Further, one often has vector quantities expressed in 
earth coordinates (such as wind and gravitational attraction) which must 
be properly introduced into the (H) coordinates. Consequently, additional 
relationships between the (H) and (E) coordinates must be derived. In 
particular, one must know how the (H) coordinate system is oriented rela¬ 
tive to the (E) coordinates, at all times. The orientation of (E) relative to 
(I) must also be known. Although other methods exist, as the Introduc¬ 
tion indicates, the use of Euler angles appears to be the most straightfor¬ 
ward for the present study. In this approach, one rotates a coordinate sys¬ 
tem, initially coincident with the (E) system, about selected axes so that 
after the rotations are performed this coordinate system will have the same 
orientation as the (H) system. One will then have the means for expressing 
vectors of one coordinate system in terms of the other. 

The reader is reminded that, in performing these rotations, the i H axis 
is to be coincident with the missile axis, while j^ remains parallel to the 

horizontal i £ - j p plane. 

To accomplish this end, assume the existence of an arbitrary vector, 
extending from the origin of the (E) system, to represent the missile axis. 
First rotate the (E) system about the k E axis by an angle 4> so that the 
i' axis (ip rotated) coincides with the projection of the missile axis on 

the i E —j K plane, as shown in Figure 1 (p 11). 
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Missile axis 



Horizontal plane 


Figure I 

One may conveniently express this rotation in matrix form, as shown in 
Appendix B. 

i' cos 0 sin 0 0 i 

j' - -sin 0 cos 0 0 ) (13) 


Another rotation is yet required about the positive j' axis, (j E rotated) so 
tha: i H (i' rotated) is coincident with the missile axis. Denote the magni¬ 
tude of this rotation by 0, leading to the matrix: 

ijj cos 0 0 -sin 0 i' 

r H ■= o i o r (i4) 

k H sin 0 0 cos 0 k’ 

To avoid difficulties later on resulting from ambiguity in the determination 
of initial conditions and in other coordinate system transformations, let 0 
be defined from 0° <_0 < 360° and -90° <_0 <.90°, for all time! Obviously 

special attention must be given to the equations of motion when 0 passes 
through in/2 since sharp discontinuities in 0 will result. Each rotation 
matrix is a linear transformation; hence, one may obtain the (H) coordinates 
directly in terms of the (E) coordinates by combining both rotations. This 
is equivalent to the following matrix multiplication. 


‘As an example. Figure 1 would require 6 to take on a negative value. 
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( 15 ) 


- 


- 



“ 


~ 


~ 


~ 

'H 


cos 

e 

0 

-sin 6 


cos p 

si n p 

0 


‘E 

4 


0 


1 

0 


-sin p 

cos p 

0 


i e 



sin 

6 

0 

cos 6 


0 

0 

1 


k E 


producing 


*H 

|h 

k H 


cos 6 cos p 
-sin tp 
sin Q cos ip 


cos 6 sin \p -sin 6 
cos p 0 

sin 6 sin p cos 6 


*E 

>E 


( 16 ) 


Note from Figure 1 (p 11) that j H (j' rotated) remains in a horizontal 
plane, or is perpendicular to k E . This is reaffirmed in Equation 16 where 
the component of j H on k E is seen to be zero. 


One may obtain the inverse of Equation 16 to obtain (E) in terms of 
(H) as follows: 


■ 


- 


• 

1 E 


cecp -sp seep 


4 

i e 


ce sp cp se sp 


47 

k l. 


-so o ce 


- k H 


( 17 ) 


For brevity cos A and sin A have been replaced by CA and SA, respec¬ 
tively. It should be noted that Equations 16 and 17 actually represent three 
scalar equations. To illustrate this, let us assume the velocity components 
of the missile expressed in the (H) system are known and we wish to ex¬ 
press the velocity in the (E) system. We have 

^ " V XH*H * V YH>H 4 V ZH k H " V XE‘k * V YE*E * V ZE k E (18 ^ 


First, obtain i„, j„, k H in terms of i E , j E , k E , which can be done by using 
Equation 16. Thus, 

i >( - C0C0i E + CB S0j i; - S0k E 

r n -S0.7; . c0r E * ok;. d 9 ) 

kj, soc^Tj ♦ * aik~ 


12 


It then remains to substitute these equations into Equation 18 and equate 
the coefficients of i E , j E , k E to V XE , V Y e. V ZE> respectively. Performing 
these operations will yield: 


V XE = VxH 0 ^'/' " V YH S ^ + V ZHS0C0 
V YE = V XH C ^ + V Ylf^ + VzH 56 ^ 


V 


ZE 


-v x „s« . v ZH c« 


( 20 ) 


This is the required transformation. 

Although much has been accomplished with Equations 13 through 20, one 
must go further to obtain the rates of change of the angles 0 and 6 as the 
missile position varies in time. 

Let Z be the angular velocity of (H) relative to (E) where, in (H) coor¬ 
dinates, 

w = “XlfH + ^YH^H + k* Ztfc H ^ 


Appendix B shows Z can also be written as 

Z - 0k E + Q\' (22) 

Using the matrices already developed, one proceeds by expressing both 
k E and j* of Equation 22 in terms of the (H) coordinates and equating the 
resulting coefficients of i H , j H , k rt to the coefficients of i H , j H , k H in 
Equation 21. This will produce 

~ -0 Sd 

WyH = b (23) 

o»zH ’ 0C0 


The rates of change of the Euler angles become simply 


W XH a, ZH 

~W ~ce 


(24) 
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One may note that several angular velocities have been introduced into 
the equations of motion, namely, <3, w, and (c3x), 1° component form, this 

amounts to nine unknowns, a completely overwhelming assignment for the 
three moment equations. Fortunately, not all these quantities are inde¬ 
pendent. It can be shown (see Ref 4, for example) that if a coordinate sys¬ 
tem (H) is rotating with angular velocity dl relative to a coordinate system 
(E), and (E) is rotating with angular velocity Q relative to a coordinate 
system (I), then the angular velocity of (H) relative to (I), is simply the 
vector sum of the individual angular velocities. Thus, 

w = Z + Q (26) 


This can be generalized to more than three coordinate systems. For 
example, continuing to use the definitions given in the Table of Symbols, 
we are also at liberty to write 

(eu-j) = Q + 6= n + 5 + Q (27) 

It is now necessary to obtain explicit expressions for 3, (w T ), and (d> T ), 
each of which is necessary to the solution of the general equations of 
motion. 

To obtain components of Q in the (H) system, one need only convert k E 
into the (H) coordinates, since the k E axis was assigned to be coincident 
with the spin axis of the earth. Therefore 

Q = Qk~ = -QS0i~ + QC0k^ (28) 


The total angular velocity, (w T ), of the missile is introduced by first 
making use of the (M) coordinate system. Since the 7^ and i f| axes are to 
remain coincident at all times, the angular velocity fl of (M) relative to 
(H) is simply; 


n 


n r 

XM M 


n r 

XH H 


(29) 


By the use of Equation 27, the total angular velocity can now be written 
in component form as follows: 


(oi T ) X H 

n XH * 

W X H - QS0 


(w T ) Yh 

W YH 


(30) 


Cti _ .. 4 

QC0 


Zll 
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One may still sense the existence of too many unknowns, namely, 0, and 
the three components of to. Fortunately Equation 23 gives a relation between 
two components of <o. Eliminating < ft produces 

gj x h ~ ~ a> ZH ^ ^ (31) 

Finally, to obtain explicit expressions for (d> T ) as given in the basic mo¬ 
ment equation, Equation 31 is substituted into Equation 30 and the latter 
is differentiated with respect to time. This leads to: 

(<VXH = ^XH- 6 i ZH tan 0 - 6 i ZH seC ^ 6 J YH- QC0 ^YH 

= 0J YH (32) 

^T^Z.H = 0J 7.H “ 2 S0W YH 

This completes the discussion. The equations needed to obtain the 
actual trajectory are summarized beginning on page 51. The next section 
considers the forces and moments which constitute the left” side of the 
equations of motion (Equations 8 and 12). 

Forces and Moments 

This section considers the forces and moments acting on the missile. 
Excluding guidance, these forces and moments carv be divided into three 
categories as follows: 


1. Aerodynamic 

2. Gravitational 

3. Jet 

Since the forces and moments are essentially as presented in Reference 
1, only a brief account of each will be given here. 


Aerodynamic Forces 


Three aerodynamic forces are considered. A brief table giving magni 
tude and direction of each is given on page 16. 
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Fore* 
Axial drag 

Normal force 
Magnus force 


Magnitude' 

P*<V' k DA 
PW ) 2 k N 

pd , (a , T) Xf /V f )k 


Di faction 

Along missile axis 

Perpendicular to missile axis 
in plane of yaw 

Perpendicular to plane of yaw J 


*(V ) is the magnitude of the missile velocity relative to the air. 

2 The yaw plane is the plane that contains both the missile axis and 
resultant velocity vector. 


These forces are illustrated in Figure 2. 



Figure 2 

The axial drag acts along the negative i|| axis, which by definition is 
directed from the CG towards the nose of the missile. Explicitly. 

(Axial drag) Xj| /hF (Y f )' k |M (M) 

The normal force components act opposite to the directions of (V ) yM 
and (V ). . Schematically we have 


1 (> 







Knowing the scalar magnitude of the resultant normal force and using 
the geometry of Figure 3, one can deduce the component of this force act¬ 
ing along j H , namely 


(Normal force) YH 


-pd J (V r ) J k N cos a 


-pd’(Vr ) 1 k N (V r ) yH 

V(v r y + (v r y 

v 'yh Z H 


(34) 


Similarly, in the k t! direction 


-pd 3 (V r ) J k N (V r ) ZH 


(Normal Force) 7 ,, 


(35) 



The Magnus force acts in a direction perpendicular to the plane of 
yaw, as indicated in Figure 4. It should be noted that, for reversed spin, 
the Magnus force will act in the opposite direction. 



Using the scalar magnitude and the geometry of Figure 4 above, it is 
readily ascertained that ir. the j H direction 


~pd 5 ((/.>_) Y ..(V )k P (V ) 7u 

(Magnus force) Y} | - ~pd 3 (<o T ) XH (V )k f siri a -I_±ll— i —_— i —— (36) 


nWyh 3 < v »>zn 


Similarly, in the k„ direction 


(Magnus force) 


pd 1 ( ( u T ) X ||(V )k F (V) Y n 


V (v y ( (v y 

' rill (711 


(37) 


Aerodynamic Munivnfi 

Proceeding in a similar vein as for the aerodynamic forces, the following 
table gives the magnitude and direction of each of the aerodynamic mo¬ 
ments considered: 


ZH 
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Moment 

Magnitude 

Direction 

Restoring 

pd J (V) J k N (r c - A n ) 

Perpendicular to the 
plane of yaw 

M agnu s 

P d> (" T )xH (V r )k F (r C " A M> 

In plane of yaw 

Yaw damping 

pd 4 (V) \/(w/ YH + (<J T )i, K k H 

Perpendi cular to 
missile axis 

Roll damping 1 

pd 4 (V f )k^ p (^ T ) XH 

Along missile axis 

Roll 2 


Along missile axis 


while for fin-stabilized missiles k^ is 

to be changed to k^. 

2 

The roll moment is used only for missiles with canted fins. 


‘V 


s used for spin-stabilized missiles, 


In obtaining the components of the above moments, one may use the 
relation "distance * force - moment,’ ’ applied in a manner consistent 
with the geometry of the system. When expressed in vector form the restor¬ 
ing moment equation becomes 

(Restoring moment) (r c - A H * [(normal force) YH j H + (normal force) ZH k H l (38) 
This produces, in component form 


(r r -AN)pd 2 (V r ) J lc N (Vr) Z H 
(Restoring momentJyjj — ^ » ' — 

v'<V)‘ v „ * (V)' z „ 


(59) 


in the in direction, and 


(Restoring moment) ZJ j 




(40) 


'< v ,»‘yh • < v .>; 


7M 


in the k,, direction. Note that both magnitude and direction are consistent 
with the above table. 

Doing likewise for the Magnus moment produces along the j|{ axis 


(Magnus moment) Ytl 


“( r C ~ A M )^»d 1 (oj T ) x v r v r )y H 


(41) 
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while along the k H axis 


(Magnus moment) 


<*r 


A M )pd 3 (a | T )x H (V )k F (V ) ZH 
Vn)y H ; ^>ZH 


(42) 


The yaw damping moment acts to diminish the yaw of the missile; 
hence, this moment acts opposite to the lateral angular velocity of the 
missile. The components of this moment are readily ascertained to be 

(Yaw damping) Y ^ - -pd 4 (V) co t ) yh ^ ^T^ZH ^ Y (43) 

(Yaw damping) ZH = ~pd 4 (V r ) \/(w T ) YH + (w T )z H k H sin 8 y (44) 

where 5 Y is shown in Figure 5 below. 



Resultant lateral 
angular velocity 

\Z(«U T ) YH + (wt^ZH 


Writing cos S Y in terms of U> T ) YH and (w Y ) ZH 


Cos S Y _ (< "T>YH 

v("J Y ) Y h 1 (wj)zu 


(45) 


and substituting into Equation 43 leads to the component along i,. 


(Yaw damping) yH ~pd*(V) (w T ) YH k (1 (46) 

similarly, for the component along FJ ( , 

(Yaw damping), () -pd 4 (V ) (w T ) 7H lt (| (47) 
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The roll damping moment acts to reduce the spin; hence, it is introduced 
into the equations of motion with a minus sign prefixed. 


The direction of the roll moment depends upon the orientation of the 
canted fins. If the cant produces positive spin (i.e., a clockwise rotation 
of the missile looking from the rear of the rocket), then this term is intro¬ 
duced with a plus sign. Like the roll damping moment, the direction is 
along the longitudinal axis of the missile. 

It should be noted that, since none of the aerodynamic coefficients 
(k DA , k N , k F , etc.) are assumed to be linear in nature, they are not to be 

taken as slopes to be multiplied by angle of attack. Rather, these coeffi¬ 
cients are point values obtained directly from experiment as functions 
both of Mach number and angle of attack. 

Table 1 (p 22) summarizes in component form the aerodynamic forces 
and moments considered thus far. 

Gravitational Fore* 

In evaluating this force, distinction must be made between spherical 
and flat earth cases. In the latter case, the gravitational force acts in a 
constant direction (independent of the missile position) while in the for¬ 
mer case, the gravitational attraction is directed towards the earth’s cen¬ 
ter from the missile CG. In both cases, specific expressions can be de¬ 
rived; the simpler case will be treated first. 

Flat Earth Ca*o 

In the absence of the earth’s rotation, 1 both theT, and k E axes are 

taken as extending in positive vertical direction. The origins of the (I) 
and the (E) coordinate systems are both located on the surface of the 
’’flat” earth. Hence, mgy (y being a unit vector representing the direction 
of the gravitational force acting on the missile) becomes. 

mgy -mgkj -mgkj.- 


tine should have little desire to reline the model to include the earth s rotation and 
yet allow the assumption of a plane to represent the earth s jjeometr>. 
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Replacing kg by its representation in (H) coordinates, similar to Equa¬ 
tion 28, one obtains 

mgy - - mgCOk^ (49) 

in a form suitable for substitution into the general equations of motion. 

Spherical Earth Case 

For this case, the direction of the gravitational force depends upon the 
missile position. In particular, let the current missile position be repre¬ 
sented by a vector R, where 

R = R Xe‘E + R YE>E + R ZE k E ^ 


as shown in Figure 6. 



Dividing R by its magnitude produces the unit vector 

R A- R XE‘K 1 R YK'E * R 7.E k k (51) 

I ft I V R xg * R yf * r ze 

Clearly, 

mgy -mgR (52) 
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Combining Equations 51 and 52 and obtaining i £ , j £ , k £ in terms of the 
(H) coordinates, there results: 


, i 


(Gravity force)xH = -- [R + R ygC0S0 - R z 1 

vR XE + r ye + R ZE 


(Gravity force) yh = -1212-[-Rxe^0 + RyeC0] (53) 

^ R XE + R YE + R ZE 

(Gravity forceJ^H = -——-[ R XeS0C0 + RyE^S0 + R ze£^1 

^ R XE + R YE + R ZE 


One should note that the gravitational acceleration, g, is a function of alti¬ 
tude. In both cases, one may write 


g 


R earth 

h 


(54) 


where 


g o - acceleration of gravity at sea level, 
R earth = ra< l‘ us °f c hc earth, 

I - current altitude of missile 
h < - (R ZE ^ or ^ at earf h case) 


+ R YE + R ZE for the spherical earth case. 


Since the gravitational force acts at the CG of the missile, there are no 
moments associated with this force. Table 2 (p 25) summarizes the results 
of this derivation in component form. 


J«t Force* 


A single jet force is considered, the thrust, which imparts forward motion 
to the missile. Ideally, this force should act along the missile’s longitudi¬ 
nal axis; however, due to imperfections in the rocke. design, the actual 
nature of propellants and other factors, there may be a component of thrust 
perpendicular to the longitudinal axis. Since the thrust vector is defined at 
a given time relative to the missile body, the thrust components will first 
be .cified in the (M) coordinate system (i.e., the coordinate system that 
is rigidly attached to the missile). 
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Using the symbols presented in Figure 7, one can easily ascertain these 
components to be as in Equation 55. 



Figure 7 


(Thrust)^ = - T cos 

(Thrust)yM = ^YM T sin 5y cos <5 A (55) 

(Thrust)^ - ^7M - T sin sin 5 A 

It is necessary, of course, to determine the components of thrust in the 
(H) coordinate system along whose axes the forces are summed. 

Figure 8 shows a typical relationship between the (H) and (M) co¬ 
ordinates. 




Figure 8 
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Here S is defined as the angle between the j M and j H axes. It is readily 
seen that 


hi 


1 0 

0 cos S 
0 -sin S 


- 


- 

0 



sin S 



cos S 


k H . 


and, from Equation 56, 


«M = *H 

j m = cos S )H + sin Sk H 
= -sin S j^j + cos S kjj 


(56) 


(57) 


Combining Equations 55 and 57 produces the required result: 


T 

1 XH 

T cos 


t yh = 

T sin cos (5^ + S) 

(58) 

t zh = 

T sin 8 j sin (5 ^ + S) 



Jet Moment* 

Two jet moments exist. One unintentionally arises from the fact that the 
line of thrust may not pass through the centroid of the missile. The other 
is an intentional jet torque which causes the rocket to spin about its axis 
of revolution. To be explicit about the former, one may define a vector? 
from the CG of the rocket to the point of application of the thrust at the 
nozzle exit. As before, this vector is defined in the missile coordinate sys 
tern, i.e., 

* r XM'M ' r YM»M * r 7.M k M 

The corresponding jet moment is now given by 

(Thrust moment) ■ r v (Thrust) 
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which in component form becomes: 


(Thrust moment) 

r YM T ZM “ r ZM T YM 

(Thrust moment) ym 

r ZM T XM ~ r XM T ZM 

(Thrust moment) 

r XM T YM “ r YM TxM 


Converting these components into the (H) coordinates and substituting 
Equation 55 finally produces 

(Thrust moment) XH - T sin ^flryM s * n ^A ~ r ZM cos ^A^ 

(Thrust momentiyn T cos ^ T ~ r XM s * n ^ T s > n 8 cos S 

- T [r sin cos 8^ - r YM cos sin S (60) 

(Thrust moment >ZH T t r ZM cos - r XM sin 6.^. sin S A 1 sin S 

\ T Ir^^ sin 8j cos 8^ - r YM cos Sjl cos S 

The final moment, the jet torque, will be assumed proportional to the 
thrust. If this moment produces positive spin, this term is to be intro¬ 
duced into the equations of motion with a plus sign prefixed. Mathemati¬ 
cally, 

(Jet torque) k T i (61) 

This concludes the discussion of forces and moments. 

One might note that V , the velocity of the missile relative to the air, 
has not been defined in the presence of winds. In the next section of this 
report explicit expressions will be derived for this vector in the (H) coor¬ 
dinate system. 

Introduction of Wind* 

This section derives expressions that can be used to introduce wind 
into the equations of motion and thus determine (V f ). For the spherical 
earth case, it will be assumed that wind data exists at preselected points 
on the earth’s surface. For a flat earth, ii is assumed that wind is pre¬ 
sented as a function of range and altitude. 

In both cases, if one knows the velocity of the missile (V M ) and the 
and the air velocity (V A ), each relative to the earth's surface, then the 
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velocity of the missile relative to the air (V ) is given by 



<Va> 



(62) 


For the flat non-rotating earth, wind is likely to be given relative to the ground; 
therefore, use of the (E) coordinate system is appropriate and Equation 62 may 
be written as 


(V ) (V x , -(V^)x|.)i{. ' (V YK -(V^)Yp)ii.- * (V/p ( V a)/.K) k f. (63) 
where (V^) is the actual wind vector. 


It is an easy matter to convert Equation 63 into the (H) coordinates. 
Doing so produces 


(V ) XH I V X1 , - (V w ) xi , 1 coco * lv YI; -(V W ) YK 1(:OS0 - |v ZK -<v w ) ZE ls/> 

(V ) Y u “(Vj-p - (\'^’)xp * l^yp (^w) Y v) f (64) 

(V,)/h l V x |.- - (V^)xE + lV Y p; - (V \/j) Y p]S0Si^ t [V Z p - (V^) Z p)C() 

Of course, (V«-) X K* (V»-)yk. and (V«')zk can be given as functions of 
range. 

The spherical earth case is a bit more complicated because of the 
geometry. Also, one must account for the motion of the air induced b_* the 
earth’s mtation in addition to wind (i.e., disturbance of the air within 
this rotating air mass). 

Recalling that Equation *i 


d.,R 


. w • H 


dt 


V Xll'll 1 v yh'u ' v /II* II 


is the velocity of the missile relative to inertial coordinates (by definition 
of u), one must also express the motion of the air induced b) the earth s 
rotation relative to inertial coordinates. By so doing one can obtain the 
....i,. c...; ;i.. ... i.... »s«. '■o'atinEair mass. Th is will be the 

equivalent of Equation 62 referred to inertial coordinates. 


The air mass rotates as a whole with the earth, this is shown physically 
by noticing that the rotation of the air mass is not apparent to one on the 
earth's surface. The velocity of this rotating air mass at the current mis¬ 
sile position (specified by the vector R) is simplv 





The velocity of the missile relative to the rotating air mass therefore 
becomes 


v M/Rot. air mass 


v - (3 * ft) 


( 66 ) 


One must now account for the wind itself. For the spherical earth case, 
winds may be conveniently introduced into the equations of motion by use 
of the azimuthal coordinate system (W). This coordinate system is de¬ 
scribed on page 6 and is illustrated in Figure 9 (p30), where the origin of (W) 
has been translated to the earth’s surface at O' for visual purposes. The 

A 

wind can be considered to be a function of the height h, of the missile CG 
above the earth’s surface, (see Figure 9). One method of designating the 
wind in the (W) coordinates is to specify wind magnitude and wind direction 
measured from North at given altitudes. It will be assumed in this report that 
the wind vector, V\(h), has no vertical component. Thus at a given altitude 
Vw(fi) shou.d be expressible only in the iw~plane, as presented in 
Equation 67, using the geometry of Figure 9. 



V) 




y' “W* 


(V w ) - |v^)| sin <£ j y + |V^ (£)l cos 4> k\f/ = (V^)y^j*jp +(^w^ZW^W (67) 

The presence of ft emphasizes the dependence of wind on altitude. The 
point O' on the earth’s surface is defined by a and y where 


a 



y -- tan 


f R 7.E 


( 68 ) 


where if 


r xe 0 


R XE 0 


r yk >0 

0 

<a <90° 

r ye <0 

270° 

_ a < 360 c 

R Y K '■ 0 

90° 

< a < 180 

r yk 0 

180° 

< a v 270 


and 


-90 ” _ y _90° 


It remains to obtain the wind components i.i the (H) coordinate system 
and add the results to Equation 66. This is done by first obtaining j w 
and in terms cf (E) and then (by the now overworked Equation 17) 
obtaining (E) in terms of (H). 

Rotating lirst about k { . through the angle « and then rotating about j' 
through the angle ft {ft -v in Fig 9). one can see that the resultant co¬ 
ordinate system will be coincident with (W). 


Performing the indicated operations results in 



CftCn C.ftSn -$ft 


* K 

-Si <i (■a 0 



Sft( it Cft 


\ . 


31 


(69) 





Combining Equation 69 with Equation 17 completes the transformation. 
Thus, 


l w 


C/SCa C/SSa -Sj8 


" 

COCi/j -S 0 S0C0 



j w 

= 

-Sa Ca 0 


C0S0 C0 S0S0 



k w 


S/3Ca S/SS a C/3 


-S(9 0 C<9 


k H 


leading to 


r 1 




“ 

*w 


^11 ^12 A,j 


5 

i 


A. A . A,, 


t 

J w 


“2 1 “2 2 “2 J 


'H 



^Jl ^12 


1 

u 1 


where 

A n = C/3CflC0C0 + CpSaCdSi/, + S/3S0 

A 13 -CfiCaSi/j + C/3S aCifi 

A ls • C0CaS0Cdf f C.ftSaSOSi/j - S/3C0 

A,, =■ -SaCOCi/j t CaC0S0 

Aj j SaSi/» 4 CaC0 

A j5 ShS9Cu . CaSDSi.) 

A„ S/JCciC0C0 t SfiSaCO S0 - C/3S.9 

A„ -SfiCaSiji • SfiSaC^j 

A„ SflCuSOCO * SflSaSflStf . C/K9 

Clearing up the remaining threads, we have 

(^ w ) (V w ) ' ( V v )^ V W ^XH*H * f V vi Jyn'n ' ^ V W^/U lc H 
where the usual algebraic manipulation results in the necessary equation 


( \ 7 XII » Aji ’ (V wV.w A >> 
( V V .. ( V V * A » ’ ( Vzw A » 
( V/.H ( Vy* A "’ ( V/.. A « 


(72) 





(73) 


Finally, the sought after vector is 

(V r ) V-(QvR)-(V w ) 

* 

Note again that Vv prescribes the wind velocity relative to the rotating 

air mass. Thus, for a missile moving on or above the surface of the earth 
with a velocity of Q ^ $ (stationary orbit or fixed surface point), in the ab¬ 
sence of winds would cause the relative velocity, V f to be zero. Writing 
Equation 73 in component form yields 


( VxH 

^XH (V\h')xh 


(Vyh 

v VH aC^Rx,, + Q^R 7h -(V w ) yh 

(74) 

l V r >ZII 

V 7.H f nS ^ R y H - ^VzH 



;here, of course, the magnitude of (V r ) is 


W V<y>* xfi 

This completes the analysis of winds. 

Initial Condition* 

Initial conditions for the spherical earth case will be derived first. 
Initial conditions for the flat earth are immediately obtainable upon spe¬ 
cializing certain of the parameters. 

In order to start the trajectory, a compete set of initial conditions must 
be provided. These are tabulated below for convenience. 


xt: 

V XH 

0(0) 

**XH 

v l- 

V Y». 

0(0) 

<|J YH 


V/H 


<■'/!! 


Those conditions that depend upon a coordinate system for representa- 
non can bo f iv« in rerms of «(.« oordinarc sysrems. ,1 des.red, as Ion* 
as known transformation equations exist between such coordinate systems. 

Of the conditions stated above, the mtssil. or.en.ation expressed in 
terms of )',<(» and »(#> is not usual!) known, per se. But the follow,n* 
measurable quantities ate (or al Ic" r 'hould bet known: 


Longitude of the launch point 
Latitude of the launch point 


A° 


B° 

Azimuthal heading of the missile G° (This will be defined more explicitly 

later in this report.) 

Angle of declination of the missile 
from the local vertical H° 

It is possible, as will be shown, to express 0(0) and 0(0) as a function 
of the variables A°, B°, G°, and H°. To do this, one rotates a coordinate 
system initially (t 0) coincident with (E) through the above measured 
angles. Upon completion of the rotations, the resultant coordinate system 
is identical in orientation with the (H) coordinate system at launch. The 
(H) system is also defined by simply specifying 0(0) and 9(0). Each of 
the resultant rotations can be expressed in matrix form, one matrix contain¬ 
ing the four known variables and the other matrix containing the two un¬ 
known variables. One is then at liberty to equate corresponding terms of 
the two matrices, since each actually represents the same coordinate sys¬ 
tem. A particular notational scheme is used to distinguish the several 
rotations performed. Attached to each unit vector is a series of primes, the 
number of which denotes the number of rotations already performed. 

The purpose of the first two rotations is to locate the i axis so that, 
when extended, it passes through the CG of the missile at time t - 0. 

This is done by first rotating about kthrough the angle A° so that the 

i' axis lies on the same longitude as the missile CG and, secondly, by 
rotating about j' through the angle B", the latitude of the missile position. 
These rotations are u^ed in the same order to obtain the position of the 
azimuthal coordinates; hence, one may use Equation 68 with « and (3 
replaced by A and B, respectively. Thus 



The purpoi-c of the next two rotations is to orient the i’ 1 axis so that it 
becomes coincident with the missile axis. 

In particular, the third rotation is abou t i' 1 through the angle G", so 
that after this rotation is performed -k’'' coincides with the projection of 


-a 


the missile axis on the j* *' — k’ 1 * plane. In matrix form, we have 



The fourth rotation is about j' 1 ' through the angle H° so that i''' ' is 
finally coincident with the missile axis. Upon examining the geometry of 
these rotations, 11° is seen to be the angle of declination the missile makes 
with the local vertical i''. 

Mathematically, 



Combining all rotations performed results in 



whe re 


b„ 

CHCBCA 

- SHSGSA - 

SHCGSBCA 

b„ 

CMCBSA 

t SHSGCA - 

SHCGSBSA 

b„ 

-C MSB - 

SHCGCB 


b,. 

-cgsa t 

SC.SBCA 


bj a 

CXiCA . : 

SGSBSA 


bj i 

SGCB 



b»i 

SIK Me A 

. CIISGSA ■ 

( 1K X.SIU A 

l*„ 

SM( MSA 

( llSt.t A • 

( IK C.SMSA 

b„ 

-SMSH • 

( IK C.( M 



Although one may feel that a sufficient number of rotations have already 
been performed to obtain i* * * ' coincident with the missile axis, there is 



yet no assurance that j"" will lie in a horizontal plane, as required. This 
situation is remedied by performing a fifth rotation about i "" through an 
angle L such that the resultant position of j"" lies in the horizontal 
plane. In matrix form 


r 


r- 



*H 


1 o o 


'fill 

in 

= 

0 CL SL 


j MM 



0 -SL CL 


k"" 


The magnitude of L is as yet unknown. Substituting Equation 78 into 
Equation 79 produces the matrix 



b„ 

bj ,CL 4 bj,SL 
-bj,SL t b s ,CL 


bn 

b 2 j CL 4 b 32 SL 
-b 2 2 SL f b 12 CL 


b 2 jCL 4 bj,SL 
-b 2 jSL 4 b,jCL 



(80) 


which enables one to convert from (H) to(E) coordinates. Equation 16 


*H 


COCO COSO -SO 


1 E 



-so co o 


»E 

k n. 


S0C0 S0S0 co_ 


k F . 


also relates the (II) and (E) coordinates. Since these sets of equations 
are equivalent at the start of the trajectory, one simply compares corres 
ponding terms, from which one obtains L as well as 0(0) and 0(0). 

In particular, 


Sm 0(0' -b,, 


Sin 0(0) -_li_ 

Cos 0(0) 


Cos 0(0) -U- 

Cos 0(0) 


(81) 


(82) 


My 



and by using 


b 23 cos L + b n sin L ■= 0 
bjj cos L i bjj sin L = C0 


(83) 


one obtains for L 


Cos L = b » 

bj j b J2 - b SJ b 3 j 

(84) 

Sin L = b :a 

b 2 a b J2 - b n b 2 j 

Although L need not be explicitly used in determining 0(0), 0(0), L com¬ 
pletely describes the missile orientation at launch. 

Since 0 was restricted between ±90°, the quadrants of each angle are 
uniquely determined in the preceding equations. 

The initial linear and angular velocities may be given in terms of either 
the (E) or the (H) coordinates, although w is usually most easily expressed 
in the (H) coordinate system. Perhaps it would be best at this point to 
illustrate the theory by a brief example. 

Let us suppose that it is desired to launch a rocket as pictured in Fig¬ 
ure 10 (p 38). We make the following assignments: 

A = 30° (Longitude) 

B - -40° (Latitude. Note the negative rotation about j' for the Northern Hemi¬ 
sphere; hence, the minus sign). 

G = 10° (190° — 180° so that the negative k'" axis will coincide with the 
line 0* P). 

H -- 35° (Angle of declination from the local vertical). 

0(0) and 0(0) can then be determined by use of the equations up to and 
including 84. 

Let us now introduce a 5° yaw angle in a plane containing the local 
vertical, measured towards the north, as pictured in Figure 11 (p 39). From 
the geometry of Figure 11 (b), one can obtain the components of the mis¬ 
sile velocity in the (E) system and use the transformation Equation 17 to 
obtain V , V , and V . fi is simply the spin of the missile as 

All iM 
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Local vertical (i " ) 








Diagram Illustrating Initial Conditions (II) 


North 


Local vertical 


Initial velocity vector 
Missile axis 



Local vertical 


Velocity vector 
Missile axis 


Figure 11 

Y) 






observed in the CH) coordinate system. and w ZH are the lateral 
angular velocities of the missile as observed from the surface of the earth. 
Sketching the rotations in sequence before numerical values are assigned 
to these quantities can be helpful. 

We will now specialize the results so that they are appropriate to the 
flat earth case. As has been mentioned previously, the vertical is to be 
the k E axis, while the ip axis will be taken as downrange. One conveni¬ 
ent way of ascertaining initial conditions is to first obtain the i" axis 
coincident with the k E axis; then, with G r and H° known, compute yi(0) 
and 0(0). Referring to Figure 10 (p 38). we replace the 30° with 0° and 
the -40° with -90° (or 270°). The i" axis thus achieves its first desired 
orientation. Mathematically, 

A 0° 

B 270° 

Our matrix Equation 78 for arbitrary angles G and H now becomes 


[" i"" 


SHCG 

SHSG 

CH 


‘E 




-sc. 

CG 

0 


! e 

(85) 

k~ 


-CHCG 

•OiSG 

SH 


. k E - 



When the appropriate quantities are substituted imo Equation 83, L be¬ 
comes zero. Using these results to obtain <£{0) and 0(0) for the flat earth 
case gives the following: 


u(0) G 

0 ( 0 ) H - 90 (86) 

L O'* 

for 

A O' 

b :?o 

The orher quantities V X(l .'‘V.H arc » as before, treated analogously, 

although R x R Y p, and R / , may (if appropriate) be assigned values of 

zero. 
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This concludes the section on initial conditions. More information will 
be provided on ' is subject when the alternate set of equations, as derived 
in the next sec. on, is used. 

Singularity Conditions ana an Altornato Sot of Equations 

This section presents an alternate set of equations to be used when 
singularity conditions are approached and, equally important, the means to 
convert to this set during flight simulation. 

As indicated under the General Equations of Motion, special attention must 
be given to the equations of motion when ft approaches * rr/2. As had been antic¬ 
ipated, the term <,y ..tan ft in the equations of morion approaches infinity as ft ap¬ 
proaches 90 . Furthermore, 0 becomes indeterminate when this singular condi¬ 
tion is reached. Unfortunately, ft 90 occurs whenever the missile axis be¬ 
comes vertical, the occurrence of which cannot be ignored. To avoid this condi¬ 
tion, the (H) coordinate system will be replaced by a new fixed-plane sys¬ 
tem labelled the (VI coordinates. As outlined in the beginning of the pro¬ 
cedure, | y is constrained to tie perpendicular to the i fe - axis (i.e., to lie in 

a specified fixed vertical plane) while iy is to be coincident with the mis¬ 
sile axis. This change necessitates redefining the Euler angles, and deriv¬ 
ing equations oi transformation between the two coordinate systems. 

To satisfy the conditions imposed upon the (V) coordinate system, the 
rotations indicated below are made. 

1. A rotation about ij.. of magnitude 4'. 

2. A rotation about j' of magnitude so that «' ' becomes the 
missile axis. 

To avoid ambiguity, the following angle restrictions were observed: 

0 V>0 

0 -ft 180 

In matrix form 



■U 




After multiplying, we have 


*v 


C0 

ses'p 

-SGC'P 


‘E 


>v 

k v 


0 

S0 

-cosy 

S Y 

C0CW 


i e 
. k E 

(88) 


Note that, as required, i v has no component along i E . 


Rather than introduce new notation, a prime over a given "omega” will 
denote that the particular angular velocity that was previously referred to 
in (H) coordinates will now refer to the (V) coordinates. 


Continuing in a manner analogous to the horizontal coordinate system, 
we can write 

«• = ’PT^ + 0“=’i’[C0l^ + S0k^] + 0j^ (89) 


or iii component form 


and for later purposes 


“kv * *» 


CO 


YV 


= 0 


(O 


ZV 


= ¥S0 


CO 


CO 


ZV 


XV 


tan 0 


(90) 


(91) 


Note that, when 0 = 0°, a singular condition is again present; however, 
investigating the geometry of the rotations for V and 0 will show that this 
occurs when the missile axis is horizontal rather than vertical, a situation 
for which the previous set of equations is applicable. 

The angular velocity of the earth in the new coordinates becomes 

Si = QiT* = [-ns0c^i“ + ns^j” - ncoc^ (92) 
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Letting ft' be the angular velocity of the (M) coordinates relative to 
(V) coordinates, we can write 


(wip - f!' + «' d' + w' +6 


(93) 


or, in component form (after using Equations 91 and 92), 


(<o'j 


ft' 


rxv xv T tan © 


zv -qsocv 


(w’) - w' + QSV 

TYV YV 


(94) 


KW’ "kv * ac0c ' 1 ' 


Differentiating Equation 94 yields 


) « Q< ^ZV 0 + “zv _ Q [C0CTO - S0SVV] 

T X V XV sin 2 0 tan 0 

(<4) yv = + ncvv 

(o^zv - <^zv - Q t s0c ’ y ® + C0SV vi 


(95) 


One may note that Equations 8 and 12, the basic component equations 
of motion (as well as most forces and moments), were derived without 
first specifying the orientation of the (H) coordinate system; hence, they 
are equally valid in the (V) coordinates. Therefore, using the "omegas” 
and the "omega-dots” derived above, along with the new transformation 
equations, produces an alternate set of equations valid for vertical orien 
tations of the missile. 

Since both sets of equations have their own singularity conditions, it 
it necessary not only to have both sets of equations available but it 
should also be possible to convert from one system to the other as the 
need arises. It is the purpose of the next few equations to establish this 
conversion. 

By definition the T axes of the (H) and (V) coordinates are coincident; 
hence, if an additional rotation is made about the i y axis so th:< j y 
rotated lies in a horizontal plane (as determined by i^. and j^. ), then the 
(V) coordinate system will become coincident with the (H) coordinates. 
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Let us denote this additional rotation by 0, whose magnitude is yet un¬ 
known. In matrix form 

10 0 i Y 

j~ = 0 CD SD j7 (96) 

k^ 0 -S<1) CD k v 

Combining the inverses of Equations 88 and 96 eventually produces 


C0 

S0SD 

S0CD 

S0SD 

CDCD - C0SDSD 

-C'VSD - 

-S0CD 

S^CD fC0CDSD 

-S*PSD 4 


Comparing this with Equation 17 yields explicit expressions for ob¬ 
taining V, 0, and D in terms o ( 6 and 0 , as follows: 

Cos 0 cos 0 cos 0 (98a 

Sin D - cos (9 sin jj . cos ip sin 6 (98b 

sin 0 sin 0 


Sin D - ~ siri V 
sin 0 


cos D ^ sin 6 cos 
sin 0 


Again no ambiguity results from the restrictions made earlier on 0 (i.e., 
the quadrants of 0, V, and D are uniquely determined). 

Knowing S', 0, and D and, of course, 0, 0, 0, and (), we are now in a 
position to determine the angular rates cf change, H* and 0. 

This is accomplished by differentiating Equations 98a and 98b, producing 

<3 SflC0 <b + CAS00 

S0 


^ sw| -sos^o * oc001 - a)S0c00 

CV sin' 0 





if ¥ happens to be either 90° or 270°, one can differentiate the second 
equation (98b), obtaining instead of Equation 100 


^ = -s0c<9<9 + sflcee 

SW sin 1 © 


(10 1 ) 


to eliminate the singularity. 

Once these derivatives are known it follows, from Equation 90, that 


o>2 V = V sin 0 


6 )' =0 
YV 


( 102 ) 


To obtain O xv , one need only equate the two expressions for (w T ) x , 
obtaining 


^ . - a S0C¥ = flxH - <u ZH tan Q - Q SS ( 103) 
tan 0 


or, since S0CW ^ S0, 


fi XV " n XH~ w ZH tan ^' 


to 


ZV 


tan 


0 


(104) 


To complete the transformation, Equation 96 again comes into use. We 
have for the invariant vectors (force, moment, and velocity) 


F' 

XV 

- F 
r XH 


L XV 

= l xh 

F' 

YV 

= F C4> - F 
YH r ZH 

S«t> 

4v 

= l yh c<I> “ l zh s4> 

FL 

= F S<t> + F„, 

C<1> 

L’ 

= L„, S$ + L , U C<1> 

ZV 

YH ZH 


ZV 

YH ZH 


V' 

- V 




XV 

XH 




V' 

- V C4> 

- V 

S<U 


YV 

YH 

ZH 



V' 

= V S<t> 

4 V 

CO) 


ZV 

YH 

ZH 



(105) 
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Other equations required for using the fixed vertical plane are 



a;, - -Sac© f Casesv 

A|, S/}CaCfl . S£S«S0S>f' - CftSQC'V 


( 108 ) 


This result and the other two components are tabulated in the Summary of 
Equations (p 51). Needless to say, (V ) is determined from Equation 72 

written in terms of the (V) coordinates. Wind for the flat earth case is 
analogous to that of the (H) coordinates. The resulting equations are also 
tabulated in the Summary of Equations section. 


Before leaving this section, three derivations are yet required. These 
involve expressing the thrust and thrust moments, specifying new initial 
conditions in the (V) coordinates, and converting from (V) back to (H). 
We begin with the thrust equations. 


Since $ is defined as the angle between the j and j axis, Figure 8 

V H 

(p 26) may be modified as follows: 


4c 



from which it is clear chat the thrust components expressed in the new co¬ 
ordinates are given by Equation 58 with S replaced by S'. Thus, 

T xv = T cos S T 

Ty V = T sin 8 tCos(§ a + S') (109) 

T zv = T sin § T sin (5 A + S') 

At changeover (i.e., when converting from (H) to (V) coordinates), angle 
S' is defined by 


S' = <t> + S 


( 110 ) 


after which the following is used 

S’ = f* O'dt + S' I 

t changeover It changeover 


Use of these equations implies that there is no reorientation of the (M) 
coordinates; hence, all quantities previously referred to these coordinates 
(such as S A , Sp r XM’ r YM> r ZM^ retna i n constant when changing from (H) 

to (V) coordinates. 

One similarly replaces S by S' in Equation 62 to obtain the thrust mo¬ 
ments in the (V) coordinates. 

This completes the transformation equations pertinent to the thrust terms. 
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Initial conditions must be provided in the (V) coordinates, because of 
the possibility of "vertical launch." Actually, much that has already been 
derived in the section on initial conditions can be used here. The only 
real change is that angle L has to be redefined relative to the (V) coordi¬ 
nates. Angles A, B, G, and H are defined as before. Let L' denote this 
new angle. It is determined by noting that j v is to have no component 

along i E , in place of kp. This is effected by equating Equation 80 (with 
H’s replaced by V’s) with Equation 88 (instead of with Equation 16 as 
was done before). Mathematically, we have 

Cos 0(0) - b u 

Sf C¥ -JjiL (112) 

s?T S0 


We have for L' 


CL* 

SL' 


b„ CV 

bj i b,j - bj, bj } 

b 3 ,CW 

b: i bjj — b,, bj, 


(113) 


where b lt , ..., b,, are defined as before. 

For the flat earth case, again setting A 0°, B 270°, we have 


Cos 0(0) SHCG 

Cos V(0) -Q1 (114) 

S0 

Sin M'(0) - -§G 

S0 


Finally, it is necessary to determine expressions for converting from 
(V) to (H) coordinates. To obtain these expressions, one can proceed in a 
manner similar to previous derivations. 

An angle f is defined such that when a rotation of this magnitude is 
performed about i v , the j v axis will then lie in the horizontal plane (i.e., 
parallel to the equatorial plane). This rotation expressed in matrix form 
will yield Equation 96, with $ replaced by <t>' . 
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Postmultiplying the inverse of Equation 88 by the inverse of Equa 
cion 96 produces 


1 


- 


~ 

‘K 


c:0 s0s<t>‘ s(-k <h’ 


*H 

•e 


30S4» efeth’ - CfisVsfl*' -CH'Sth’ - COSl'C^' 


ill 

k E 


-SfiC'f' SH'C <ti' , C0CYS<D' -S4>' SV + C0C<fc'( V 


. k H. 


( 115 ) 


Comparing Equation 115 with Equation 17 below 



coco 

-Sy) 

soco 

easy) 

co 

soso 

-so 

0 

CO 



results in the following transformations: 


SO S0C¥ 


CO 

C0 

sc- 

S0S¥ 

CO 

co 

s<h' 

-SO 

c<t>' 

soco 


S0 


S0 


( 17 ) 


(116) 


No ambiguity results, since -1 _o •. 


To obtain the angular velocities, one differentiates Equation 116, 
obtaining 


-) CHCVH--SHSW 

CO 

CO cos' 0 


If 0 is 90 or 770 ’, one can use 


u >s00 - c.Qsoo 

Si,'/ cos' 1 0 


( 117 ) 


( 118 ) 


( 119 ) 


Equation 104 completes the transformation of the angular velocities. 


40 




Solving for QxH^ eac ^ s to 


fiXH Oxv + w 7H tan ^ t ——— (120) 

tan 0 

where both oy M and are defined by Equation 23- 

When Equation 96 is used, the forces, moments, and velocities can be ex 
pressed in (H) coordinates, given their representations in (V) coordinates. 
These auxiliary equations are tabulated in the following section of this 
report. 

This, therefore, completes the transformation from the (V) to the (H) 
system. A complete tabulation is presented in the following section. 

Summary of Equations 

The equations of motion as derived in the preceding five sections are 
summarized below. In the writing of the equations, the primes have been 
dropped when working in the (V) system. No ambiguity should result. 


SO 



Force Equotions (H-Coordinotes) 
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tan 6 - 8S0I ly w y(l - rL» Y(J I x “ <'> z „ tan ^ “ 8S011 




where 
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II. For The Flot Earth Cose 
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Equations of Transformation 
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SHCBSA - CHSGCA + CHCGSBSA 





CHSB - SHCGCB 


/ 
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= 0 if trajectory starts in H-Coordinates, T = T c hang eo ver otherwise. 




In all equations 
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Fore* Equations (V-Coordinates) 
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0 if trajectory starts in V-Coordinates, T = ^changeover 




Equations of Tronsfofmotion (V) -* (H) 
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In all equations 




RESULTS AND DISCUSSION 


The final equations resulting from the derivation are tabulated in the 
preceding section. There are, however, many aspects of the equations that 
have not been discussed. It is the purpose of this section to consider 
briefly some of these aspects and limitations and to mention what is in¬ 
volved in incorporating certain refinements in the equations as presently 
derived. 

One might begin by noting that a typical rocket trajectory consists of 
several flight phases. Explicitly, these phases, typical for a two-stage 
rocket, may be tabulated as follows: 

1. Acceleration of booster and main stage 

2. Coasting of booster and main stage 

3. Separation of booster and main stage 

4. Coasting of main stage 

5. Acceleration of main stage 

6. Free flight of main stage. 

These phases are in direct correspondence with those used in the six- 
degree-of-freedom digital computer trajectory program at Picatinny Arsenal. 
The user is thus able to run any or all of these phases. Phase 3 is effected 
by imparting a short thrust to the main stage, simulating the actual booster 
separation. 

Since the aerodynamic coefficients, ccnte«s of gravity, pressure, etc. are 
not necessarily identical for all phases, several sets of such data are 
often required to simulate a complete trajectory. The equations in the pre¬ 
ceding section of thL report a.e used tor ail phases. 

To numerically solve the differential equations of motion, a modification 
of the third order Runge-Kutta technique is used (Runge-Kur >Gill). This 
numerical scheme appears to be adequate for both low - and high-spin pro¬ 
jectiles, although inclusion of thrust malalignments remains to bt investi¬ 
gated. As is indicated in the introduction, however, i: appeared that th ; s 
numerical scheme was not satisfactory when high angular 'ates of change 
were implicit in the equations themselves. This condition pertains to high- 
spin projectiles for which the fixed-plane coordinate system is rigidly 
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attached to the missile (as was originally done). This is an advantage of 
the present coordinate system. 

Other aspects of these equations, such as the use of thrust modification 
factors, the atmospheric model (ARDC Atmosphere of 1959). etc., will be 
discussed in a forthcoming report. 

Some of the refinements of the equations are discussed below. 

1. Introduction of an ellipsoidal earth model. This modification has 
far-reaching effects on the equations. For example, the direction and mag¬ 
nitude of the gravitational attraction (consistent with the earth model) 
would vary with missile position. Additional coordinates would have to be 
specified on the surface of the ellipsoid to properly introduce wind data. 

In addition, an iterative scheme is necessary to compute the altitude of the 
missile if it is naturally defined as the shortest distance from the missile 
CG to the earth model. Initial conditions would also be modified accordingly. 
For increased ranges and accuracy requirements of trajectory simulation, 
the ellipsoidal earth refinement may well become necessary. 

2. Treatment of asymmetrical missiles. Asymmetrical missiles can 
arise from two sources. The more severe situation is present if the ex¬ 
ternal missile configuration does not possess rotational symmetry. For 
this case, the transverse moments of inertia could vary about the fixed- 
plane coordinates as the missile rotates. Also, specific data relating the 
resultant point of application of the aerodynamic forces would have to be 
provided for accurate results. A second asymmetrical condition can occur 
for externally symmetrical missiles when the missile CG is offset by a 
prescribed amount from the longitudinal axis of the missile. This obviously 
modifies the lever arms in the moment equations and may also introduce 
time-dependent moments and products of inertia. 

3. Other refinements could include guidance and launcher effects, both 
of which can play prominent roles in trajectory analysis. 

In addition to the limitations listed earlier, two assumptions were im¬ 
plicit in the derivation. The first is that the motion of the earth about the 
sun was neglected; for most sub-orbital trajectories this phenomenon can 
be ignored. The second assumption is that the rate of change of inertias 
was neglected (i.e., i terms during the thrusting periods). For missiles 
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possessing excessively high burning rates with large angular velocity com¬ 
ponents, the magnitude of these terms should be investigated. 

A final word about the equations is that there is no estimate of the dis¬ 
persion of the missile. This requires the computation of several trajectories, 
each for a slightly different initial condition, with appropriate statistical 
combinations of the various ranges and deflections from a fixed standard. 
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APPENDIX A 


Rotating coordinate systems 

It is the purpose of this Appendix to express the rate of change of a 
vector relative to fixed coordinates, in terms of rotating coordinates. Briefly, 
this analysis answers the following question: "Given a particle whose 
motion is known, how would the particle’s motion appear to someone 
situated on a coordinate system that is itself in rotation?” 

Obviously, one must know the motion of the rotating coordinates relative 
to the fixed coordinates. To be explicit, let (I) denote fixed coordinates 
and (H) a rotating coordinate system, whose origins are both coincident. 
Further, let <ibe the angular velocity of (H) relative to (I). This motion 
might be as depicted in Figure A-l, where the three dotted discs represent 
the paths of each axis of the (H) system for the given constant u. 



Fig A-l 
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Let Rbe an arbitrary vector. Clearly, R has a representation in the (H) 
system as follows: 


R = R XH ‘H + R YH >H + R ZH k H 


(A-l) 


To determine the rate of change of ft relative to the (I) coordinates in terms 
of the (H) coordinates one must account not only for the changing magnitude 
of ft, but also for the variation of the unit vectors i H , j H , k H (which help 

represent ft) relative to (I). Mathematically: 


UiA , , 

— = r xh ‘h + r yh >H + R ZH k H + r xh ~r 
dt dt 


dt 


dt 

(A-2) 


The first three terms of Equation A-2 define simply how ft itself is changing 
independent of any moving coordinates. The latter three terms describe the 
motion of (H) with respect to (I). It is noteworthy that the I subscript de¬ 
notes differentiation with respect to the (I) coordinates. To elaborate fur¬ 
ther, if I were replaced by H then Equation A-2 would read 



dt 


R XH ‘h 


r yh ’h 


R ZH k H 


(A-3) 


since, for purposes of the differentiation, i H , j H , k H would be fixed. 

• ▼ d lc 

The task remains of obtaining expressions for . . . , ““jp* To do 

this, let us "isolate” the i vector as presented in Figure A-l, and define 
some new quantities as given in Figure A-2. Further, let 
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r. . 


i H (t + At) denotes the position of i H at time t + At. Clearly, by the defi¬ 
nition of a derivative, 


i H (t + At) - i H (t) 


Liin . . 

dt At-»o At 


and from the geometry of Figure A-2 

i H (t + At) - i H (t) = A i H ~ i H sin 0 [«At] 

we may write 


(A-4) 


(A-5) 


d I ‘H 


i^ sin 0 «At 

dt 


At 


and in the limit 


jljjj 

dt 


= Lim 
At-0 


i H sin 0 oAt 

At 


i H sin 0 M 


(A-6) 


Equation A-6, however, is very reminiscent of the cross-product, namely 

(A-7) 


d i r H 


dt 


= MX 1 


l h 


and, if one examines the directions, one can indeed state 
d l*H - T 


dt 


*mx i H 


(A-8) 


Since i H is typically representative of each coordinate axis, we also have 


*H - r 

— a M X 


di ku ■* ~ 

JL” = MX L 

dt ^ 


(A-9) 
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It is natural to substitute Equations A-8 and A-9 into Equation A-2. This 
produces 


djR • — • — * — -» — 

~ 7 ~ = r xh ‘h + r yh *h + r zh ^ + r xh ( MX 


+ Ryu ( wx i H )+ Rzh (“ x k H^ 


which one may write as 


^—= ^4— + wx r xh ‘h + “ x r yh Jh + “* Rzh k H 

dt dt 

or finally 

d.R d„R - -* 

—!— = —— + u x R 

dt dt 

which, in the text proper, is the basis for Equation 4. 


(A-10) 


(A-11) 


(A -12) 
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APPENDIX B 

Matrix representation of rotations 

Assume the geometry of Figure B-l, where a rotation ofmagnitude 0 has 
been performed about the k £ axis, producing new vectors i' and j’. 



Fig B-l 


Clearly the projection of i' on the i E axis is i' cos 0 = cos 0 

(all vectors shown are of unit magnitude). Continuing in like manner for 

all possible combinations, one may form the following table: 

‘i >1 

i* cos 0 sin 0 0 

j* -sin 0 cos 0 0 

k' 0 0 1 

Here, for example, -sin 0 is to be interpreted as being the projection of j' 
on the ie axis. 

If we now choose to represent i' in terms of its components along i , j , 
and k£, we write 

i* * cos 0 i E ♦ sin 0 jg ■» 0 kg (B- 1) 
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Though matrix multiplication, we can write not only i', but also j* and 
k', not separately as in Equation B-l, but together in one matrix similar 
to the given table. Thus; 


i' 


cos 0 sin 0 0 


’ »E 

i' 

= 

- sin 0 cos 0 0 


*E 

k' 


0 0 1 


.y 


or Equation 13 of the text. 

It might be mentioned that this development is not intended to be rigorous; 
however, it should enable one to obtain these rotation matrices without 
difficulty. 

Finally, let us complicate the above and now assume that i' and j* are 
rotating about kg with some angular velocity u,. Clearly u, is directed 
along the kg. axis, and further, in time At we have the relation 


i«,i ■ 



(B-3) 


or, in the limit, as At approaches zero 


!w,l 


dt 


(B-4) 


Finally, combining both magnitude and direction, we can write 

«• (B*5) 

Ve can do likewise for a rotation about j'. obtaining, for example. 


«*, «i' (B-6) 

where 9 is an angle defined analogously in the i* - k* plane as 0 was de¬ 
fined in the i e - jg plane. 

Since angular velocities may be added, we can combine Equations B-3 
and B-6 and write 

ii*' • 

tit in | ♦ hii 0 kg • ti j' (B-7) 

which, in the text, is Equation 22. 
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